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Abstract 

The concentration of primary radiation induced defects has been previously esti- 
mated considering both the explicit mechanisms of the primary interaction between 
the incoming particle and the nuclei of the semiconductor lattice, and the recoil 
energy partition between ionisation and displacements, in the frame of the Lind- 
hard theory. The primary displacement defects are vacancies and interstitials, that 
are essentially unstable in silicon. They interact via migration, recombination, an- 
nihilation or produce other defects. In the present work, the time evolution of the 
concentration of defects induced by pions in medium and high resistivity silicon 
for detectors is modelled, after irradiation. In some approximations, the differen- 
tial equations representing the time evolution processes could be decoupled. The 
theoretical equations so obtained are solved analytically in some particular cases, 
with one free parameter, for a wide range of particle fluences and/or for a wide 
energy range of the incident particles, for different temperatures; the corresponding 
stationary solutions are also presented. 
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1 Introduction 



A point defect in a crystal is an entity that causes an interruption in the lattice 
periodicity. In this paper, the terminology and definitions in agreement with 
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M. Lannoo and J. Bourgoin [1] are used in relation to defects. 

Vacancies and interstitials are produced in materials exposed to radiation in 
equal quantities; they are the primary radiation defects, being produced either 
by the incoming particle, or as a consequence of the subsequent collisions 
produced by the primary recoil 

In previous papers [2,3] we calculated the concentration of radiation induced 
primary defects (CPD) in semiconductors exposed to hadron irradiation. Due 
to the important weight of annealing processes, as well as to their very short 
time scale, CPD is not a measurable physical quantity. 

In silicon, vacancies and interstitials are essentially unstable and interact via 
migration, recombination, annihilation or produce other defects. The system 
evolves toward natural equilibrium. 

The problem of the annealing of radiation induced defects in semiconductor 
materials is not new. Several models, empirical or more theoretic, have been 
previously proposed to explain these phenomena. Damask and Dienes devel- 
oped, in the 1960s, a theoretical model that considers, as distinct cases, the 
kinetics of vacancy-interstitial annihilation with interstitial trapping at im- 
purities [4], of vacancy-interstitial annihilation with interstitial trapping at 
impurities [5], and of vacancy-interstitial annihilation with diinterstitial for- 
mation [6]. In the first two cases, analytical solutions have been obtained by 
the authors, while in the last case only numerical solutions are possible. The 
escape of interstitials and vacancies from clusters and their subsequent interac- 
tions have been modelled by Tsveybak et al. [7]. Phenomenological approaches 
of the annealing could be found in Moll et al. [8] and Bates and co-workers [9] 
and references cited therein. 

Most of the old studies were dedicated to electronic devices, made from low 
resistivity (in the Qcm range) silicon. Recently, the behaviour of the silicon ma- 
terial with medium and high resistivity (tens flcm - kVtcm) in intense hadron 
fields started to present very much interest, these materials representing a 
major option for the detectors at the new generation of particle colliders and 
for space applications. 

In the present paper, the time evolution of the concentration of primary defects 
induced by pion irradiation in silicon is studied in the frame of a simplified 
phenomenological model based on direct interactions between the primary 
induced defects and the impurities present in the material. 
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2 Equations for the kinetics of radiation induced defects: general 
formulation 

The starting hypothesis of the physical model is that the irradiation process 
is unique at the determined particle fiuence, and that the time of irradiation 
is shorter that all characteristic times of annealing. 

In the frame of the model, equal concentrations of vacancies and interstitials 
are supposed to be produced by irradiation, in much greater concentrations 
than the corresponding thermal equilibrium values, characteristic to each tem- 
perature. Both the pre-existing defects and those produced by irradiation, as 
well as the impurities, are assumed to be randomly distributed in the solid. 
An important part of the vacancies and interstitials annihilate. The sample 
contains certain concentrations of impurities which can trap interstitials and 
vacancies respectively, and form stable defects. 

In the present paper, vacancy-interstitial annihilation, interstitial migration to 
sinks, vacancy and interstitial impurity complex formation as well as divacancy 
formation are considered. The sample could contain more impurities that trap 
vacancies or interstitials, and in this case all processes are to be taken into 
account. 

The following notations are used: V - monovacancy concentration; / - free 
interstitial concentration, J\ - total impurity " 1" concentration (impurity " 1" 
traps interstitials and forms the complex C\)\ C\ - interstitial-impurity con- 
centration: one interstitial trapped for one complex formed; J 2 - total impu- 
rity "2" concentration (impurity "2" traps vacancies and forms the complex 
C2); C2 - vacancy- impurity concentration: one vacancy trapped for one com- 
plex formed; V 2 - divacancy concentration. All concentrations are expressed 
as atomic fractions. 

The chemical scheme for reaction rates has been used, and the differential 
equations corresponding to the above picture are: 



dV 
~dt 



K X VI - K 5 V ( J 20 - C 2 ) + K 6 C 2 - K 7 V 2 + K 8 V 2 



(1) 



dl 
~dl 



K X VI - K 2 I - K 3 I ( J 10 - d) + K 4 d 
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i ^ = K i I(J 1Q -Ci)-K A Ci (3) 
at 



^ = K 5 V (J 20 - C 2 ) - K 6 C 2 (4) 



The bimolecular recombination law of interstitials and vacancies is supposed 
to be a valid approximation for the present discussion, because at the con- 
centrations of vacancies of interest, only a small fraction of defects anneal by 
correlated annihilation if their distribution is random (see the discussion in 
[10]). 

If N is the total defect concentration, expressed as atomic fraction: 



N = V + / + 2V 2 + Ci + C 2 (6) 
then it satisfies the differential equation: 



dN 

f - -mVI (7) 



The initial conditions, at t — 0, are: at the end of the irradiation, there are 
equal concentrations of interstitials and vacancies: Iq — Vq] the concentrations 
of impurities are Jio and J 2 o respectively, there are no complexes in the sample: 
Cio — C 20 = 0, and no divacancies V 20 = 0. 

The reaction constants K\ and K 3 are determined by the diffusion coefficient 
for the interstitial atom to a substitutional trap, and therefore K\ = K 3 : 



K x = 30z/exp {-En/k B T) (8) 
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where En is the activation energy of interstitial migration and v the vibra- 
tional frequency. The reaction constant K 2 (related to interstitial migration 
to sinks) is proportional to the sink concentration a: 

K 2 = av\ 2 exp {-E iX jk B T) (9) 
with A the jump distance. 

K4 characterises the decomposition of a complex into an impurity and an 
interstitial, and is given by: 

K 5 and K 6 describe the formation and decomposition of vacancy-impurity 
complexes, while K 7 and K 8 the formation and decomposition of divacancies. 
and are given by: 

K 5 = 30vexp(-E i2 /k B T) (11) 
with E i2 the activation energy for vacancy migration; 

k ^<- e W) (12) 

K 7 = 30uexp{-E i2 /k B T) (13) 



K > = ^(-n£r) (14 » 

where B\ is the binding energy of C±, B 2 the binding energy of C 2 and B 3 is 
the corresponding binding energy of divacancies. 
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Due to the mathematical difficulties to solve analytically the complete differ- 
ential equation system, some simplifications are necessary. 



3 Hypothesis, approximations and discussions 

The interstitials are much more mobile in silicon in respect to vacancies, and 
are characterised by an activation energy of migration a factor of two times 
smaller. 

This difference between the activation energies En and E i2 respectively justi- 
fies the introduction of two time scales, and the separate study of the processes 
involving interstitials and vacancies respectively. In the first stage, vacancy 
interstitial annihilation and interstitials migration to sinks are studied. The 
concentration of interstitials decays much rapidly than the concentration of 
vacancies. A cut-off condition for / is imposed (a p times decrease of inter- 
stitial concentration, p being an adjustable parameter of the model). The 
vacancy concentration determined by this procedure is the initial value for 
the processes comprised in the second stage, and will be denoted by Vmuiai- 
The formation of complexes by vacancies is considered less important, and is 
neglected in the following discussion. 

So, for the "first stage", after some simple manipulations, from (1) and (2), 
the equation: 




K 2 



(15) 



has been obtained, with the solution: 




(16) 



and: 



* = 




(17) 



K 2 ln(V(t)) 
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Imposing the cut-off condition for the concentration of interstitials, both V in i tia i 
and the characteristic time could be found. As specified, Vinitiai is used as initial 
vacancy concentration for the second step in the analysis, where the equations 
(4) and (5) are considered. Vmuiai depends on Vb and p, and it is temperature 
independent. This system of equations, expressing the kinetics of divacancy 
and vacancy-impurity formation, have no analytical solution. These processes 
are governed by the initial concentrations of vacancies and impurities. 

If the impurity that traps vacancies is phosphorus, the limiting cases cor- 
respond to low initial doping concentration (high resistivity, uncompensated 
materials) and very high impurity concentration (low resistivity), respectively. 
In both cases the equations could be decoupled, and analytical solutions are 
possible. 

So, if the formation of vacancy- impurity complexes is not so important, the 
main process responsible for the decay of vacancy concentration is divacancy 
production. In this case, the time evolution of the vacancy concentration is 
described by the equation: 



v M = m\-* +R TT ^<J(i \ (18) 

where: 



R = y/K 8 (K 8 + 8K 7 Vi) (19) 
while the increase of the divacancy concentration is given by: 



The stationary solution for V(t) is given by: 



For n-type high doped Si, the process described by eq. (4) is the most probable. 

7 



(20) 



If Jo is the initial concentration of impurities, and the initial concentration of 
complexes is zero, than: 



with: 



R* = yjKi + Kl (Vi - Jo) 2 + 2K 5 K 6 (ft + J ) (23) 
and with the stationary solution: 




and 



c = v; - y (25) 



4 Results and physical interpretations 

In Figure 1, the absolute value of the CPD per unit fluence (CPD) induced by 
pions in silicon is presented as a function of the kinetic energy of the particle. 
This curve has been obtained in the hypothesis listed in references [11,12]. In 
the kinetic energy range between 100 to 300 MeV, pion - nucleus interaction 
is resonant in all waves, the A 33 resonance is produced and a pronounced 
maxima in the CPD is observed at approximately 150 MeV followed by a slight 
monotonic decrease at higher energies. Some local maxima are also presents 
at high energies, but with less importance. 

The process of partitioning the energy of the recoil nuclei (produced due the 
interaction of the incident particle with the nucleus, placed in its lattice site) 
in new interaction processes, between electrons (ionisation) and atomic motion 
(displacements) has been considered in the frame of the Lindhard theory. 
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The CPD multiplied by the fluence is the initial value of the concentration of 
vacancies and interstitials, and in the forthcoming discussion it is expressed, 
as specified, as atomic fraction. Without loss of generality, will shall consider 
that the primary defects are produced by 150 MeV kinetic energy pions. 

The following values of the parameters have been used: En = 0.4 eV, E i2 = 
0.8 eV, £i= 0.2 eV, B 2 = 0.2 eV, B 3 = 0.4 eV, v = 10 13 Hz, A 2 = 10 15 cm 2 , 
a = 10 10 cm" 2 , The value of the parameter^ is taken 100. 

The time evolution of the concentration of vacancies, normalised to the con- 
centration of vacancies created by irradiation is represented in Figure 2. Flu- 
ences in the range 10 11 — 10 15 pions/cm 2 have been considered. The weight of 
the annihilation process in respect to interstitial migration to sinks increases 
abruptly with the fluence. At low fluences, up to 10 11 pions/cm 2 , the annihila- 
tion has a low importance. With a good approximation, it could be considered 
that after 0.2 sec. after irradiation the vacancy concentration is saturated. As 
could be seen from the figure, observable effects in V/Vq appear for fluences 
higher than 10 12 pions/cm 2 . Experimentally, for silicon particle detectors irra- 
diated with hadrons, some modifications in the electrical characteristics have 
been observed in the fluence range 10 12 — 10 13 cm -2 [13,14]. 

All the curves correspond to room temperature (20 °C) annealing. 

The interval of time after which the concentration of interstitials decreases 
p times, given by eq. (17) with V(t) = V ini ti a i is a characteristic time of the 
processes in the first stage, and depends on temperature and fluence. 

In Figure 3, this characteristic time is represented as a function of fluence, in 
the fluence range 10 11 — 10 15 pions/cm 2 , for three temperatures of interest for 
high energy physics applications: -20, and +20 °C It may be seen that, up 
to 10 13 pions/cm 2 , the characteristic time is independent on the fluence, and 
for higher fluences a decrease of the characteristic time is to be noted and a 
fluence dependence is obtained. 

The initial concentration of vacancies for the second stage is V initia i previously 
calculated. The second stage starts at a time equal to the characteristic time. 

Two limiting cases have been studied, low doping concentration (high resis- 
tivity uncompensated material), and low resistivity, phosphorus doped silicon. 
The formation and decomposition of divacancies is a good approximation for 
the first, while the second is related to V — P complexes. 

The time evolution of the concentration of vacancies when divacancy forma- 
tion is considered is illustrated in Figure 4a. As underlined in Section 3, the 
time scale of this process is much longer in respect to the first stage. The 
origin of time in Figure 4a is the beginning of the second stage, this means a 
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characteristic time after the start of the annealing. 

The curves in this figure are normalised to the concentration of vacancies 
remained after the first stage. Values of the initial concentrations of vacancies 
that correspond to fluences in the same range as in the previous analysis have 
been considered. It could be observed that the time after which V(t)/V in i tial 
saturates is in the range 10 7 — 10 8 sec, and increases with the irradiation 
fluence. All curves have been calculated also for 20 °C temperature. 

The time evolution of divacancy concentration is represented in Figure 4b, 
also in a similar manner with the vacancy concentration in Figure 4a. 

The case of high initial impurity (phosphorus) concentration, where divacancy 
formation is neglected in respect to complex formation, is presented in Figures 
5a and 5b respectively. 

Three initial concentrations for phosphorus in silicon (expressed as atomic 
fractions) have been considered: 10~ 5 , 10~ 7 , and 5-10~ 10 , corresponding to 0.1 
flcm, 10 Qcm, and 1 kflcm resistivity respectively. The time evolution of the 
concentration of vacancies is represented in Figure 5a, while the corresponding 
one for the concentration of complexes in Figure 5b, both normalised to the 
concentration of vacancies remained after the first stage. An irradiation with 
10 13 pions/cm 2 has been considered. 

The origin of time in Figures 5a and 5b is displaced in respect with the origin 
of the annealing, in the same way as in Figs. 4a and 4b. 

The curves corresponding to 10 11 — 10 15 cur 2 are superimposed on the same 
curve, with the exception of the tail at times longer than 10 4 sec. where the 
values corresponding to higher fluences are slightly higher. 

The time interval after which V(t)/V init iai saturates is shorter than in the case 
of divacancy, and at room temperature is in the range 10 3 — 10 4 sec. for the 
doping concentrations considered. It shortens with the increase of impurity 
concentration. 

Some explicit considerations must be done about the formation of vacancy- 
impurity complexes. The mechanisms supposed above can be used both for 
boron and for phosphorus impurities. In this case, the corresponding processes 
are [8]: 



Bi + V -> [Bi - V] (26) 

and respectively: 
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P s + V - [P 9 - V] 



(27) 



While the complex formed by boron is unstable and self anneals bellow room 
temperature, the interaction between a V and a P s leads to the formation of 
an E centre which is stable in the same conditions [7]. 

Interactions between interstitial oxygen (another very studied impurity in the 
last time) and free vacancies is described as a higher order process (third order 
in the [8] and fourth power in [15]). If the process is stopped as a first order 
one, the time evolution of the concentrations is not different from the case of 
phosphorus, studied before. 

If two or more impurities that trap vacancies are considered as existing si- 
multaneously in silicon, the system of coupled equations must be solved. Only 
numerical solutions for particular cases are possible. 

If two or more impurities that trap vacancies are considered as existing si- 
multaneously in silicon, the system of coupled equations must be solved. Only 
numerical solutions for particular cases are possible. 

So, in the frame of this model, three characteristic times have been found 
for 20 °C self annealing in silicon: one in the order of 0.2 sec, corresponding 
to vacancy - interstitial annihilation and interstitial migration to sinks, one 
around 10 4 sec, related to vacancy-phosphorous formation and decomposition, 
and another one around 10 7 sec, divacancy related. 

The present model could be applied to the annealing of primary radiation 
defects induced in Si by any other hadron, the difference coming from the 
value of the CPD. 



5 Summary 

The time evolution of the primary concentration of defects induced by pions 
after irradiation process has been modelled. 

In this model, vacancy-interstitial annihilation, interstitial migration to sinks, 
vacancy and interstitial impurity complex formation as well as divacancy for- 
mation have been considered. In the considered hypothesis, the time evolution 
of impurity concentrations has been decoupled into two steps, the first one in- 
volving vacancy-interstitial annihilation and interstitial migration to sinks, the 
second vacancy-complex and divacancy formation. 
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The equations corresponding to the first step have been solved analytically 
for a wide range of particle fluences and for a wide energy range of incident 
particles, and for three different temperatures: -20, and +20 °C. 

The approximations that permit to decouple the differential equations repre- 
senting the time evolution processes in the second step have been studied, and 
the processes have been treated separately. The simultaneous consideration of 
more processes in this last step of the model is possible only numerically. 

The predictions of the model are in general agreement with the experimental 
results, and suggest a possible connection with the self annealing mechanisms. 
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Figure captions 



Figure 1: Energy dependence of the absolute value of the concentration of 
primary defects per unit fluence induced by pions in silicon. 

Figure 2: Concentration of vacancies created by irradiation versus time, nor- 
malised to the concentration of vacancies for the first stage, for different par- 
ticle fluences, expressed in cm-2. 

Figure 3: Characteristic time of the first stage versus the irradiation fluence, 
for -20, and +20 °C temperature. 

Figure 4: a) Time evolution of the concentration of vacancies, normalised to 
the concentration of vacancies remained after the first stage when divacancy 
formation and decomposition are considered, for different particle fluences, 
expressed in cm" 2 , b) Same as (a) for the concentration of divacancies. 

Figure 5: a) Time evolution of the concentration of vacancies, normalised to the 
concentration of vacancies remained after the first stage when V — P complex 
formation and decomposition are considered. The three curves correspond 
to different initial impurity concentrations, expressed in atomic fractions, b) 
Same as (a) for the concentration of complexes. 
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